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We describe a generic mechanism by which a system of Dirac fermions in thermal equilibrium 
acquires electric charge in an external magnetic field. To this end the fermions should have an 
additional quantum number, isospin or color, and should be subject to a second magnetic field, 
which distinguishes the isospin/color, as well as to a corresponding isospin chemical potential. The 
role of the latter can be also played by a nontrivial holonomy (Polyakov loop) along the Euclidean 
time direction. The charge is accumulated since the degeneracies of occupied lowest Landau levels 
for particles of positive isospin and anti-particles of negative isospin are different. We discuss two 
physical systems, where this phenomenon can be realized. One is monolayer graphene, where the 
isospin is associated with two valleys in the Brillouin zone and the strain-induced pseudo-magnetic 
field acts differently on charge carriers in different valleys. Another is hot QCD, for which the 
relevant non-Abelian field configurations with both nonzero chromo-magnetic field and a nontrivial 
Polyakov loop can be realized as calorons - topological solutions of Yang-Mills equations at finite 
temperature. The induced electric charge on the caloron field configuration is studied numerically. 
We argue that due to the fluctuations of holonomy external magnetic field should tend to suppress 
charge fluctuations in the quark-gluon plasma and estimate the importance of this effect for off- 
central heavy-ion collisions. 

PACS numbers: 13.40.-f, 72.80.Vp, 12.38.Mh 



I. INTRODUCTION 

Interactions of chiral fermions with magnetic fields are 
the origin of a number of nontrivial physical phenomena, 
which have been intensively studied in the last few years 
[1] . Some of the most notable examples of such phenom- 
ena are the Chiral Magnetic Effect (CME) [2, 3], which 
is the generation of electric current in a system of chiral 
fermions at finite chiral chemical potential, and the Chi- 
ral Separation Effect (CSE) [4, 5], the generation of chiral 
current in a system of chiral fermions at finite chemical 
potential. The most striking feature of such phenomena 
is that electric current is generated due to the magnetic 
field, so that no electric field is needed to initiate the 
transport of charge/chirality. Rather, this current has a 
topological nature related to the axial anomaly for chiral 
fermions and spin projection in the lowest Landau level 
[2, 6] and flows without dissipation even for interacting 
fermions [5, 7, 8]. Experimentally, such phenomena are 
supposed to manifest themselves in specific features of 
the distributions of charged particles produced in heavy 
ion collisions [9]. 

In this paper we describe a similar phenomenon, in 
which an electric charge, rather than current, is gener- 
ated in a system of Dirac fermions in the presence of an 
external magnetic field. For this reason we call it "charge 
catalysis by a magnetic field" . To this end the fermions, 
in addition to charge and chirality, should carry yet an- 
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other quantum number, which we call isospin T3 (not to 
be confused with the term isospin commonly used in par- 
ticle physics) and assume that it takes values T3 = ±1. 
In QCD, this additional quantum number is naturally 
associated with color. The fermions should also inter- 
act with an isospin magnetic field, or a chromomagnetic 
field in the case of QCD, for which the vector potential 
A'^ enters the action as it/jr^A' j^tp. We consider a grand 
canonical ensemble at nonzero isospin chemical potential 
fi3 which enters the action as HzipT^Qij} an d thus favors 
particles with T3 = +1 and anti-particles with T3 = — 1 
and disfavors particles with T3 = — 1 and anti-particles 
with r 3 = +1. For QCD, the role of imaginary isospin 
chemical potential can be also played by the time-like 
component of non-Abelian gauge field. The generated 
charge is also imaginary in this case, and shows up only in 
charge fluctuations after the functional integration over 
the non-Abelian gauge field. While the mechanism of 
the described effect is not directly related to the axial 
anomaly, it is still similar to CME in the sense that charge 
generation (which also manifests itself as spatial separa- 
tion of positive and negative charges, as we will demon- 
strate in Subsection II C) is caused by magnetic, rather 
than electric field, and has a purely quantum origin re- 
lated to the two-dimensional index theorem. As we argue 
in Subsection II B, for a classical gas of charged particles 
this effect would be absent. 

It should be stressed that it is the assumption of a 
grand canonical ensemble that makes the generation of 
charge possible: the system is not isolated, but rather 
coupled to an environment, with which it can exchange 
particles. Thus the extra charge comes from the envi- 
ronment, and charge conservation is not violated. The 
system then acquires electric charge only if all three fac- 
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tors are present: an ordinary magnetic field, an isospin 
magnetic field and a nonzero isospin chemical potential 
/i 3 . At vanishing temperature, the latter can be arbitrar- 
ily small. The mechanism of this phenomenon can be 
also reversed, so that an isospin polarization is generated 
at nonzero ordinary chemical potential /x (which couples 
to electric charge) , again in the presence of both ordinary 
and isospin magnetic fields. 

We should note here that this work was in fact orig- 
inally motivated by an attempt to describe the Chiral 
Magnetic Effect [2, 3] by using calorons, topologically 
nontrivial classical solutions of finite-temperature Yang- 
Mills equations [10-12], in external magnetic fields as a 
model for topological transitions in QCD at finite tem- 
perature (see Subsection IV B). While we have found no 
traces of electric current flowing along the magnetic field, 
we have noticed the emergence of a nonzero imaginary 
charge density induced by a magnetic field parallel to 
the caloron axis. The present work was initiated as an 
attempt to understand this peculiar effect. 

The paper is organized in the following way: in Section 
II we start with a simple description of charge catalysis at 
zero temperature and at small values of the isospin chem- 
ical potential, and in Subsection II B discuss the case of 
arbitrary temperatures/chemical potentials. The depen- 
dence of the induced charge on isospin chemical potential 
turns out to be quite nontrivial and is in some sense simi- 
lar to Shubnikov-de Haas oscillations. In Subsection II C 
we consider the case of spatially localized magnetic and 
isospin magnetic fields. We demonstrate that even in the 
case when the net induced charge is zero, there is still an 
excess of charge in the region with nonzero field strength, 
which is compensated by an excess of charge of opposite 
sign in a "halo" around that region. Thus in such a 
setting the magnetic field induces spatial charge separa- 
tion. In Section III we illustrate the described mecha- 
nism on a very simple and at the same time physically 
relevant example of strained monolayer graphene. In this 
case the isospin index corresponds to two distinct valleys 
in the graphene Brillouin zone and the isospin magnetic 
field is the pseudo-magnetic field arising due to mechan- 
ical strain [13, 14]. In Section IV we discuss possible 
manifestations of charge catalysis in finite-temperature 
QCD, for which the role of isospin is played by the color 
of quarks and that of the isospin chemical potential by 
a nontrivial holonomy (Polyakov loop) of a gauge field 
configuration. Wc first consider a simple example of a 
constant non-Abelian gauge field and demonstrate that 
while the net electric charge vanishes in this case due to 
integration over all values of the holonomy, the described 
phenomenon still manifests itself in the suppression of 
charge fluctuations. In Subsection IV B we consider a 
more realistic example of an SU (2) caloron gauge field 
configuration. We show numerically that the induced 
charge is localized between the monopole and the anti- 
monopole which constitute the caloron, that is, in the 
region where the chromo-magnetic field is approximately 
constant. Finally, in Subsection IV C we roughly esti- 



mate the possible contribution of the discussed effect to 
charge fluctuations in off-central heavy-ion collisions, us- 
ing both the model of a uniform chromomagnetic field 
with fluctuating holonomy (Subsection IV A) as well as 
a model of a dilute caloron gas [15]. 



II. CHARGE CATALYSIS IN UNIFORM 
MAGNETIC FIELDS AND ISOSPIN CHEMICAL 
POTENTIALS: MECHANISM 

A. The limit of zero temperature and small isospin 
chemical potential 

Let us denote the strength of the ordinary magnetic 
field as B, and the strength of the isospin magnetic field 
as F. For simplicity we first assume that the fermions 
are (2 + l)-dimensional, but our arguments can be easily 
generalized to (3 + l)-dimensional case. The magnetic 
and the isospin magnetic fields in (2 + 1) dimensions have 
only a single component, B = B xy — d x A y — d y A x = 
const and F = F xy — d x A' y — d y A' x = const. The energy 
spectrum of fermions in such a background is given by 
the Landau levels 



= y/2\qB + T 3 qF\n, 



(1) 



where n = 0, 1,2, . . . and the labels q = ±1 and T3 = ±1 
denote particle/antiparticle states and the isospin quan- 
tum number, respectively. The degeneracies of these lev- 
els are 

g n \q$B + T 3 q$> F \ (2) 

where go — 1, g n >o = 2 and $b = (27r) 1 J ' dxdy B and 

5 

$>F = (2tt) 1 J dxdyF are the fluxes of the magnetic 

s 

and the isospin magnetic fields through the area S of 
the system. For definiteness, we have assumed periodic 
boundary conditions for spatial directions, which implies 
that $s and $f take integer values [16]. 
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FIG. 1: Schematic illustration of the degeneracies of the low- 
est Landau levels and the Fermi energies for particles and 
anti-particles with T3 = ±1 for <&p = 3, <f>B = 1. The number 
of short horizontal lines in each column denotes the degener- 
acy of the lowest Landau level for the corresponding quantum 
numbers and the long horizontal line denotes the zero energy 
level. 
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At nonzero isospin chemical potential ^3 the Fermi en- 
ergies for particles/anti-particles with T3 = ±1 become 

Ep' T3 ^ — (7x3/13. For simplicity let us first assume that 
/i 3 < min \J2\F -B\, ^2\F + B\j and that the tem- 
perature is zero (these conditions will be relaxed in the 
next subsection) . In this case only the lowest Landau lev- 
els for particles (q = +1) with t 3 = +1 and anti-particlcs 
(q = —1) with T3 = —1 will be occupied. The degen- 
eracies of these levels are |$s + 3> F | and | — &b + &f\, 
respectively and thus different (sec Fig. I). Therefore, if 
both <&b and <E> F are present, there is an unequal num- 
ber of particles and antiparticles on these levels, and the 
system has a nonzero electric charge: 

Q = \$ B +® F \-\-® B + $ F \ 

= 2sign($ F $ s )min(|$ B |,|$ F |). (3) 

If the strength of the isospin magnetic field F is kept 
constant, the charge grows linearly with B for \B\ < 
\F\ and then saturates at Q = ±2|$ F |. If \B\ < \F\, 
increasing the flux $s by one flux quantum produces two 
elementary charges. We also note that since the existence 
of zero energy level in the magnetic field is topologically 
protected in (2 + 1) dimensions, the result (3) is also valid 
for non-uniform gauge fields as long as /13 is smaller than 
the first non-zero energy level. 

Finally, we note that the mechanism described above 
can be reversed by switching on the ordinary chemical 
potential which couples to electric charge and thus in- 
creases the Fermi energy for particles and decreases it 
for anti-particles. Then only the lowest Landau lev- 
els for particles will be occupied, independently of r 3 . 
The degeneracies of the lowest Landau levels for par- 
ticles with T3 = +1 and T3 = — 1 are again given by 
\&b + &f\ and I — $b + ®f\i respectively and thus 
different again. Therefore, a nonzero isospin charge 
Q3 = \$b + $ F |-|- ( E > s + < I , F |is produced. 



B. Charge catalysis at arbitrary temperatures and 
isospin chemical potentials 

Let us now discuss how the generated charge (3) is 
changed at higher values of isospin chemical potential 
and/or at nonzero temperature, when higher Landau lev- 
els are also populated. Consider first the case of zero 
temperature but an arbitrary value of isospin chemical 
potential. The equation (3) is then easily generalized to 
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where [• • -J is the floor function. The first line coming 
from the lowest Landau levels is the same as in Eq. (3). 
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FIG. 2: Occupied Landau levels for particles in the T3 = +1 
sector (leftmost column on Fig. 1) and anti- particles in the 
T3 = — 1 sector (rightmost column on Fig. 1) at some finite 
value of an isospin chemical potential. The degeneracies per 
area are labeled below as [ ± B + F\. Horizontal dotted lines 
illustrate two values of /i3, one "small" (i.e. ^3 < 2 — B + F\) 
and one large. 



The second line in (4) counts the charge of the particles 
on higher Landau levels (left side of Fig. 2 and the left- 
most column in Fig. 1), while the last line comes from the 
contribution of anti-particles (right side of Fig. 2 and the 
rightmost column in Fig. 1). The factor of two in the sec- 
ond and the third lines comes from the spin degeneracy 
9n>o = 2 of the higher Landau levels. 

The expression (4), although being conceptually very 
simple, leads to quite a nontrivial dependence of the gen- 
erated charge on the chemical potential and the mag- 
netic field, see Fig. 3 for an illustration. Depend- 
ing on how many Landau levels occupied by particles 
or anti-particles lie below /13, the charge Q can take 
both positive and negative values. The charge decreases 
if one of the more closely spaced negatively charged 
states sinks below the Fermi energy and increases if 
this happens for one of the more widely separated posi- 
tively charged states. It is also easy to prove that the 
charge (4) is a periodic function of /i§ with a period 
equal to A/i§ = ^LCM(|$ F + $ B |, |$ F - where 
LCM (n, m) stands for the least common multiple of two 
integers n and m. 

The dependence of the induced charge density on the 
magnetic field at fixed isospin chemical potential exhibits 
oscillations similar to Shubnikov-de Haas oscillations, see 
Fig. 3 below. However, while in the latter case the period 
of oscillations tends to zero at zero magnetic field, here 
this happens at F close to B, when the spacing between 
energy levels for anti-particles with T3 = — I becomes 
very small. 

In case of finite temperature the floor functions in (4) 
are smoothed out by the corresponding Fermi factors, so 
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FIG. 3: Top: the density (a) = Q/S of generated charge as a 
function of isospin chemical potential squared (in units of the 
next-to-lowest Landau level E\ = 2|S- F\) for (F-B)/(F + 
B) — 1/6 (blue solid line, red dashed line and yellow dotted 
line) and (F - B)/(F + B) = 5/13 (green dot-dashed line) 
and different temperatures. Bottom: the density of generated 
charge at /13 — 1.5 V~B as a function of isospin magnetic field 
F for different temperatures. 



that the charge is given by 
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We have evaluated these sums numerically. The resulting 
dependence of the generated charge on /j, 3 , F and B is il- 
lustrated in Fig. 3. One can see that now the dependence 
of the generated charge on /i3 , F and B becomes contin- 
uous. However, the charge still changes sign as ^3 is 
changed. Since the asymptotic density of states for Lan- 
dau levels in both sectors with r 3 = ±1 is equal (since 
smaller spacing between Landau levels is compensated 
for by smaller level degeneracy), at very high temper- 
atures and/or isospin chemical potentials the generated 
charge should vanish. Asymptotic equality of the density 
of states in both sectors implies also that the described 



effect should vanish in the classical limit, when a macro- 
scopically large number of levels is occupied. 

C. Spatial charge separation 

As discussed in the Introduction, the generation of 
charge is only possible for the grand canonical ensem- 
ble, that is, for a system coupled to some environment, 
from which the charge is "borrowed" , and charge conser- 
vation is not violated. In this Subsection we investigate 
in more details the situation in which the charge den- 
sity is generated locally, but the net charge is still zero. 
We show that this situation can be interpreted as spatial 
charge separation. 

Let us formally assume that the space is infinite and 
the magnetic field strengths are described by some lo- 
calized functions 1 . In this case normalizable eigenmodes 
with zero energy exist also for non-integer values of the 
fluxes $s, $f [17]. Their number is given by the floor 
functions of $b and <&f, so that the equation (3) reads 
in this case 

Q= L|$f + *b|J - LI -*f + $b|J • (6) 

In deriving (6) we have again assumed that there is a 
small but finite isospin chemical potential /13 so that only 
the particle zero energy states with T3 = +1 and antipar- 
ticle zero energy states with T3 = — 1 are filled. 

We now restrict ourselves to the situation in which 
both fluxes <& p and $ b are nonzero, but the total charge 
(6) is still zero. For example, one can take §p = n + 1/2 
and < $s < 1/2. If the region in which the magnetic 
field is localized is sufficiently small, one can still have 
very large magnetic fields locally. We would like to study 
here whether the charge is generated locally in this case, 
and how it is distributed. 

As the simplest example let us assume that both fields 
B and F take constant values in the region r < R, where 
r is the radial coordinate in two spatial dimensions, and 
are equal to zero outside of this region. The wave func- 
tions of the zero energy states can be found analytically 
as [17] 

i^^^e-^Ht n = 0,l,...,JV ± , (7) 

where z = x + iy, s + is the unit spinor which is the 
eigenstate of the projection of the spin operator on 
the magnetic field with eigenvalue s = +1/2, N — 
LI^f i ®b I J — 1 and <fi(x,y) is the solution of the Poisson 
equation 

(dl + d 2 y ) <t>± (x, y) = B (x, y)±F (x, y) . (8) 

It is easy to see that the eigenvectors (7) are orthogonal 
to each other for cylindrically symmetric distributions of 



This means that they decay sufficiently quickly at infinity 
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magnetic field (in the absence of cylindrical symmetry 
this might not be the case). Solving the equation (8) for 
magnetic field which is localized in a cylinder of radius 
R we obtain 



0± (r) = 



(B ± F) (r 2 - R 2 ) /4, 
(B ± F) R 2 /2 \n{r/R) 



The wave functions (7) are then: 
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Although the normalization constants can be calcu- 
lated analytically, the resulting expressions are not very 
instructive and we do not give them here. It is im- 
portant to note, however, that for large |$b ± ®f\ and 
n <C ± <frp | they behave as 



1 

7^ 



\F±B\ 



-|* B ±*J7| 



(11) 



Under the assumption that all the zero energy states 
are occupied, the charge density is given by 



N+ 
n=0 



(12) 



where AT ± = [\& B ± $f|J - 1- This result is illustrated 
on Fig. 4, where we plot the charge density as a function 
of the radial coordinate r for the values of $b and $f 
for which the net generated charge (6) vanishes. One 
can see that while positive charge is accumulated in the 
region with nonzero field strength (r < R), the excess 
of negative charge is accumulated in the exterior of that 
region. At large distances r R, the negative charge 
density behaves as 



cr(r) ~l/r 2+2e , 

e =min (|$ s ± <&J 



|*B±MI), 



(13) 



where the last expression is just the smallest fractional 
part of the fluxes | ± $ F | . Thus the charge distribution 
has a negative "heavy tail" which compensates for the 
positive charge accumulated in the region with nonzero 
field strengths (this tail is difficult to see on Fig. 4 since 
we have only shown the central region with r ~ R). 

Let us now consider charge density for r <C R and as- 
sume that the fluxes | ± | are large. Since at small 



r the eigenfunctions (10) behave as i/v (r) 



the con- 



tributions of states with n ^ 1 are strongly suppressed, 
and one can safely extend the summation over n in (12) 
to infinity. Using now the asymptotic expression (11) for 
the normalization constants c^, we find for the charge 



FIG. 4: Charge density (12) as a function of the polar coor- 
dinate r. The fluxes are $s = 0.4 and $f = 1.5,2.5,10.5 
for solid-blue, dashed-red and dot-dashed green lines, respec- 
tively. 
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which is is precisely the charge density which can be ob- 
tained from the expression (3) valid in the case of uniform 
fields. We thus conclude that even if the system does 
not exchange electric charge with the environment, the 
charge generation still takes place due to spatial charge 
separation, and for sufficiently large fluxes |$b±<I>f| the 
charge density in the center of a region with nonzero field 
strength is given precisely by (3). One can see from Fig. 4 
that the charge density near r = is indeed very close to 
this asymptotic value (which is equal to a R 2 = 2$s/7r 
if <f>_B < $f) for $_f > 2.5. For smaller values of $f the 
charge density is somewhat higher near the origin and 
smaller at the periphery. The requirement of sufficiently 
large flux simply means that the typical extent of the 
wave functions of the Landau states Rb ~ 1/ yB should 
be much smaller than R, so that the fermions effectively 
do not feel the boundary of the region with nonzero F 
and B. 

Another interesting example is that of magnetic fields 
with the following profile 



B{r) = 



^bpI 
r 2 +P%f 
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24>fP 2 f 
(r 2 +p 2 F f 



(15) 



for which the zero energy eigenstates again take the form 
(10) with the functions ± (x,y) given by 



± (r) = * B /2 ln(l + r 2 /p|) 
±$ F /2 ln(l + r 2 /p F ) ■ 



(16) 
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FIG. 5: Top: charge densities, Eq. (17), in units of p\ as a 
function of r for <&b = 0.4, pf/pb = 0.01 and &f = 1.5 (solid 
red line), = 2.5 (dashed blue line) $j? = 3.5 (dot-dashed 
green line) and $f = 10.5 (dotted yellow line). The vertical 
thin dashed line indicates the value of pf- Bottom: absolute 
value of the same charge densities on a logarithmic scale and 
for a much larger range of r. 



Note that for unit flux this profile is the profile of 
the winding number density of the soliton in the two- 
dimensional O (3) model, where this density can also be 
interpreted as the magnetic field of an auxiliary Abelian 
gauge field [18]. 

The charge density in this case reads 
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where we have assumed that Qp > In Fig. 5 we plot 
the charge density (17) for the case of a strong F field 
localized in a small region with pp <C Pb- It is interest- 
ing to note that in this case the charge density has a dip 
in the center (although remaining positive for r < pp). 
This dip emerges due to the suppression of the densi- 
ties of eigenstates with large n by a factor proportional 
to r 2n . In the case of uniform fields this contribution 



conspired with the gaussian envelope to make the dis- 
tribution uniform. Here, however, the envelope is not 
gaussian and the effects of higher modes is to create a 
dip in the center of the region with nonzero strength of 
isospin magnetic field F . Again, one can see that the 
excess of positive charge around r = is compensated 
by negative charge density at large r, which again decays 
as (13). In order to illustrate this slow decay of negative 
charge density, on the second plot on Fig. 5 we plot the 
absolute value of a (r) at large r ^> pp on a logarithmic 
scale. One can see that indeed at some large r ^> p F the 
charge density changes sign and then decays very slowly. 



III. CHARGE CATALYSIS IN STRAINED 
GRAPHENE 

Let us demonstrate now how the described phe- 
nomenon can be realized in graphene, a two-dimensional 
hexagonal crystal lattice of carbon atoms. The disper- 
sion relation of valence electrons in graphene has two 
zeros - the Dirac points K + and K~ . The vicinities 
of these points are called valleys K + and K~ , respec- 
tively. The low-energy excitations in both valleys can be 
described as Dirac fermions [13], which couple to mag- 
netic fields in the usual way. However, by applying me- 
chanical strain in certain directions one can also induce 
a pseudo-magnetic field which has opposite signs for the 
two valleys, but does not distinguish between particles 
and holes (playing the role of anti-particles) [14]. This 
pseudo-magnetic field has been recently studied experi- 
mentally, and Landau levels which correspond to pseudo- 
magnetic field strengths larger than 300 Tesla have been 
observed [19]. 

Thus it is natural to associate the isospin quantum 
number with different valleys. Then we should also in- 
troduce the "valley" chemical potential pk, so that the 
two valleys have different numbers of charge carriers in 
them. A practical method which can be used to gener- 
ate valley currents (and thus imbalance in valley popula- 
tions in some part of the sample) by cyclic application of 
strain and chemical potential has been discussed recently 
in [20]. Note that in contrast to the general derivation 
given above, the degeneracies of the lowest Landau lev- 
els in such a setting are given by \qB + t 3 F\ and the 
chemical potential couples to T3 rather than to 5T3 (see 
Fig. 6, top panel), but the principle by which the charge 
is generated is the same. Namely, now the lowest Landau 
levels for particles and holes in the valley K + are pop- 
ulated. Again, the level degeneracies for particles and 
holes in this valley are different and given by |$s + $p| 
and I — $s + $f I ) respectively. Thus there are more parti- 
cles than holes, and the graphene sample should acquire 
nonzero electric charge. 

Considering the realistic experimental situation in 
which the strength of the magnetic field is much smaller 
than the strength of the pseudo-magnetic field [19], the 
charge density per unit area can be estimated from (3) as 
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field will then change the degeneracies of energy levels 
in different color sectors, as in Eq. (3). Recently, it was 
demonstrated that this difference of level degeneracies 
leads to the appearance of nonzero electric dipole mo- 
ment in the background of topologically nontrivial con- 
figurations of non-Abelian gauge fields [2'2). However, the 
existence of the "color chemical potential" , which would 
be required to produce electric charge in this situation, 
is clearly prohibited by gauge invariance. 

Still, the weights of different Landau levels in the parti- 
tion function can be different for different color sectors if 
the gauge field configuration in the Euclidean path inte- 
gral representation of the thermal partition function has 
a nontrivial non-Abelian holonomy (Polyakov loop) 



FIG. 6: Schematic illustration of the degeneracies of the 
lowest Landau levels and the Fermi energies for particles 
and holes in different valleys in the graphene Brillouin zone 
at nonzero "valley" chemical potential (top panel) and at 
nonzero chemical potential (bottom panel). As in Fig. 1, the 
number of short horizontal lines in each column denotes the 
degeneracy of the lowest Landau level for the corresponding 
quantum numbers and the long horizontal line denotes the 
zero energy level. 



Ql (SB) ~ 10 11 cm _2 T _1 . For sufficiently strong mag- 
netic fields of order of 1 T, such charge density can be eas- 
ily measured using modern experimental methods such as 
scanning tunnelling spectroscopy. 

Again we can invert the mechanism of charge accu- 
mulation: if instead of the valley chemical potential 
the ordinary chemical potential is nonzero, the lowest 
Landau levels will only be populated for particle states 
in both K + and K~ , as illustrated in Fig. 6 at the 
bottom. The degeneracies of these levels are given by 
|$ B + $ F | and | $b — | . Thus there are more particles 
in K + than in K~ , and the "valley charge" with density 
\Qk+ ~ Qk- 1/ (SB) ~ 10 11 cm- 2 !- 1 is generated. Esti- 
mating the value of charge density in the valley pumping 
device of [20] as \Qk+—Qk- \/S ~ 10 8 cm -2 (the ratio of 
the estimated "valley charge" per pumping cycle to a typ- 
ical area of the device) , we see that for B ~ 1 T the mag- 
netically induced valley asymmetry is much larger. Thus 
the described phenomenon can be also potentially useful 
for applications in (at present hypothetical) valley-based 
electronics [20] . We note that the possibility of breaking 
the symmetry between the two valleys due to combined 
effect of mechanical strain and magnetic field has also 
been discussed in [21], but the possibility of charge catal- 
ysis was not considered. 



IV. HOT QCD IN MAGNETIC FIELD 

In QCD, we associate the isospin quantum number 
with the color of quarks, and the isospin magnetic field 
with the chromomagnetic field. The external magnetic 



P (x) = V exp 



dx°A% (x°,x) 



(18) 



around the compact Euclidean time direction x° £ [0, (3\. 
Here V is the path-ordering operator, j3 = T^ 1 is the 
inverse temperature and r Q are the generators of the 
gauge group. The time-like component of the gauge field 
Aq (x°,x) couples to quarks in a way similar to imag- 
inary "color" chemical potential. This is easiest to see 
after diagonalizing P, which then does not couple the 
states with T3 = +1 and T3 = —1 anymore. However, 
in the Euclidean path integral Aq plays the role of the 
Lagrange multiplier for the Gauss law and thus should 
be integrated over. By virtue of charge conjugation in- 
variance the net electric charge is clearly zero after such 
integration. Nevertheless, as we demonstrate below, the 
"virtual" charge which appears due to non-trivial holon- 
omy still manifests itself as a negative contribution to the 
fluctuations of electric charge. 



A. Charge catalysis in constant chromomagnetic 
fields and the effect of holonomy 

As a model assumption, let us consider the configura- 
tion of SU (2) ® U (1) gauge fields in (3 + 1) dimensions 
with uniform and parallel magnetic and chromomagnetic 

i 7 !? and a non- 



fields with strengths B = B xy and F 
trivial holonomy with a constant 



A o (x° 



T3 



(19) 



In the language of spontaneous symmetry breaking, v 
plays the role of a Higgs vacuum expectation value. 

Such a field configuration can be thought of as an ap- 
proximation of the gauge field in the so-called caloron so- 
lution, which is the topologically nontrivial saddle point 
of Euclidean path integral of finite-temperature gauge 
theories and for gauge group SU (2) consists of a pair of 
non-Abelian monopole and anti-monopole [11, 12]. The 
chromomagnetic field between them can be roughly ap- 
proximated as constant. While direct numerical studies 



of induced charge in the caloron background will be de- 
scribed in the next Subsection IV B, here we use this 
simple approximation to illustrate the basic features of 
charge catalysis. 

Assuming periodic boundary conditions in spatial di- 
rections, the fermion partition function for massless 
quarks in such a background can be written as 



]DgZf(v,li)/V 



dk z 
~2tt 



+ 00 

E 



9n X 



log(l + e ^<^+^/ 2+ M + 
2n \ J 

B ~ Fl log (l + e -^<*.-Wa-^) + 



2tt 
\B-F\ 

2tt 
B + F 



log 1 + e 



-/3E~ , -iPv/2+Pfi 



+ I ~ I + i I log (l + e -PB- hm ^pv/^ I , (20) 

where k z is the momentum along the magnetic fields, V 
is the total spatial volume and 



E n,k z 



^2\B±F\n + k 2 z 



(21) 



are the Landau levels for particle/anti-particle states 
with r 3 = ±1. We have also introduced an infinitcs- 
imally small chemical potential /i as a source for elec- 
tric charge. In the end it will be set to zero. The 
the third and the fourth lines in (20) correspond to the 
particle and anti-particle states in the t 3 = +1 sec- 
tor, while the fifth and the sixth line are the same for 
the T3 = — 1 sector. To simplify the discussion further, 
let us assume that the temperature is sufficiently small 
(T 2 < min(|F - B\, \F + B\)), so that the lowest Lan- 
dau levels with E Q k = \k z \ dominate in (20). Calculat- 
ing now the derivative of (20) with respect to fi at fi = 
and evaluating the resulting integrals over k z analytically, 
we obtain for the expectation value of the charge density 
over the fermion states: 



T d 
(q)f ^y 2 } 1 -Q^f 



\B + F\-\B-F\) 



(22) 

f3v E [-2vr,27r] 



This result is very similar in form to (3), except for the 
factor of i in front. The origin of this factor is clear from 
(20): the holonomy plays the role of imaginary chemical 
potential for the T3 quantum number and thus induces 
imaginary electric charge. However, in contrast to (3), 
the quantity (22) does not correspond to the expectation 
value of a Hermitian operator in the ground state of some 
Hamiltonian. The correspondence with the canonical for- 
malism is only recovered after integration over the holon- 
omy, which obviously results in zero expectation value 
of the physical charge operator. Consider, however, the 



fermionic expectation value of the squared charge: 

(Q 2 ) f = T 2 ZJ 1 (v, fi) JL^Zf (v, /x) (23) 

_ VT max (F, B) Av 2 V 2 min (B 2 , F 2 ) 
~ ^ (2^ ' 

The first summand is the contribution of connected 
fermion diagram and is independent of the holonomy. 
The second term corresponds to disconnected fermion di- 
agrams and is exactly the square of the "virtual" charge 
given by (22). We see thus that the effect of holonomy 
is to reduce the charge fluctuations. We note that a sup- 
pression effect of nontrivial holonomy on charge density 
at finite chemical potential (without magnetic fields) has 
also been found in [23]. 

Finally, we have to integrate the above expression over 
all values of the holonomy with the weight Zt (v, \i = 0). 
Assuming again that the lowest Landau level gives the 
dominant contribution to the partition function and per- 
forming the integration over k z in (20), we obtain 



Zf (v, fj, = 0) = const x exp I — 



Vmax(\FL \B\ 



8vr 2 T 



(24) 



Approximating the compact integral over v G 
[— 2irT, 2irT] by a Gaussian integral over the whole real 
axis (which makes sense for sufficiently large volumes), 
we obtain the following result for the expectation value 
of the squared charge 



:q 2 ) 



VTmax(\F\, \B\ 



max OF 2 , £ 2 ) J ' ^ 1 



Thus charge fluctuations become suppressed as the mag- 
netic field grows until the strength of chromomagnetic 
and magnetic fields become equal. At that point the fluc- 
tuations become zero, and then start growing again. Note 
also that for the integration over the holonomy in (25) 
we did not use the SU (2) Haar measure dvsin 2 (w), 
but rather the U (1) Haar measure on the interval 
[— 2ttT; 2nT] . We did so since for the caloron configu- 
rations the holonomy and the chromomagnetic fields be- 
tween the caloron constituents have the same orientation 
in color space. If, however, one decides to take into ac- 
count different color orientations of chromomagnetic field 
and holonomy, the expression (20) should also be modi- 
fied, since Landau levels in different color sectors would 
no longer be the eigenstates of the operator A^T a . We do 
not consider these complications here. 



B. Electric charge and current on caloron 
configurations 

As a more realistic model of quantum fluctuations in 
non-Abelian gauge theory one can consider a gas of topo- 
logically nontrivial gauge field configurations, such as an 
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FIG. 7: Characteristic gluonic quantities of an 577(2) caloron 
with maximally nontrivial holonomy on a 24 3 x 6 lattice in two 
spatial dimensions with monopoles separated along X3. The 
third color component of the chromomagnetic field is depicted 
by arrows. The Polyakov loop is characterized by its contour 
lines 1/2 ■ TrP = ±0.4. At the monopole/anti-monopole the 
Polyakov loop becomes trivial 1/2 ■ TrP — > ±1, while out- 
side of the monopoles the Polyakov loop becomes diagonal 
(proportional to T3) and traceless, 1/2 • TrP — > 0, v — > nT. 



instanton gas [24]. However, at finite temperature in- 
stantons are no longer the solutions of classical Euclidean 
Yang-Mills equations. A generalization of the instanton 
solution to Euclidean space with compact time direction 
is the caloron solution, either with a trivial [10] or non- 
trivial [11, 12, 25-27] holonomy at spatial infinity. The 
former is stable at sufficiently high temperatures, while 
the latter is believed to dominate at lower temperatures 
[15]. Much like instantons, calorons are strongly localized 
solutions. In SU (2) , they consist of a pair of Bogomol'ny- 
Prasad-Sommerfeld monopolc and anti-monopole, sepa- 
rated by some distance. In the region of space between 
the monopole and the anti-monopole the chromomag- 
netic and chromoelectric fields can be roughly approx- 
imated as constant and parallel to the caloron axis, see 
Fig. 7. Hence, if the magnetic field is also parallel to the 
caloron axis, such solution can realize the case of paral- 
lel magnetic and chromomagnetic fields described above. 
We can thus expect that the magnetic field superimposed 
on the caloron configuration leads to the appearance of 
nonzero (imaginary) electric charge. 

Since explicit expressions for the Dirac propagator in 
the combined background of a caloron plus magnetic field 
are not known, we have to find it numerically. To this 
end, we put the SU (2) caloron field configurations on 
the lattice. We use calorons of maximally nontrivial 



holonomy meaning that the asymptotic Polyakov loop 
becomes traceless (i.e. v = ttT asymptotically), for which 
the monopole and anti-monopole have equal topological 
charges 1/2. 

Since we are studying the phenomena related to quark 
chirality, we use the chirally invariant Neuberger's over- 
lap Dirac operator [28] to calculate the vacuum expec- 
tation value of the current. Since the caloron configu- 
ration has a unit topological charge, the Dirac operator 
has exactly one zero mode, which makes it formally non- 
invertible in the massless limit. To regulate this trivial 
divergence, we introduce a small quark mass m q = 0.01 
in lattice units. Then the Dirac operator is inverted using 
the SHUMR method [29, 30]. In order to check that our 
results are not affected by finite- volume and boundary ef- 
fects, we consider caloron configurations of the same size 
on 16 3 x 4, 20 3 x 4 and 24 3 x 4 lattices. The dependence 
on the lattice spacing is checked by rescaling the caloron 
configuration that has been put on the 16 3 x 4 lattice by 
a factor 3/2 and putting it on the 24 3 x 6 lattice. Thus we 
assume that the lattice spacing is a = 1 for the 16 3 x 4, 
20 3 x 4 and 24 3 x 4 lattices, and a = 2/3 for the 24 3 x 6 
lattice. Since the current and charge densities are oper- 
ators of canonical dimension 3, their expectation values 
are also multiplied by (3/2) 3 . Similarly, the magnetic 
field strength is rescaled by a factor (3/2) 2 . Note that 
the magnetic flux is a dimensionless quantity which is 
not changed when we simultaneously increase the lattice 
size by some factor and decrease the lattice spacing by 
the same factor. Of course, there is no intrinsic physical 
dimensionful scale in our problem, since we deal with a 
fixed gauge field configuration. The physical scale can 
be restored, for example, from the characteristic size of 
calorons in the phenomenological model of the caloron 
gas [15]. This will be done at the end of the paper, where 
we will discuss possible experimental consequences of our 
results. 

To start with, we consider the imaginary part of the net 
electric charge Jo (xq) and the net electric current Jb (a;o) 
in the direction of the magnetic field B as a function 
of Euclidean lattice time xq. Spatial components of the 
current are manifestly real, since in Euclidean space they 
are still given by functional derivatives of the free energy 
over the gauge field A )1 (x): 



U (*) )/ ~ logdet (V [A„ (x)}) 



-iTr (V-^A^x)}^), i = l,2,3, (26) 



) f we again denote the expectation value 
to the fermion fields only. Since the 



where by ( . . . 
with respect 

Dirac operator V [A^ (x)] = 7 P (<9 M — iA^ (x)) is anti- 
Hermitian and Euclidean gamma-matrices are Hermi- 
tian, the resulting expression is real. On the other hand, 
the expectation value of the charge is the derivative of the 
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FIG. 8: Imaginary part of the net electric charge Jo (xq) (left) and the net electric current in the direction of the magnetic field 
Jb (xo) (right) as a function of the discrete lattice time Xq on the 16 3 x 4 lattice. The data points are slightly displaced along 
the horizontal axis in order to make the plots more illustrative. 



partition function with respect to the chemical potential: 

(Q)f ~ ^logdet(D[^(a;)] + /i7o)|^o 

~ Tr (V- 1 [A^ (x)} 70 ) , (27) 

for which the real part is identically zero at /i = 0. As dis- 
cussed above in Subsection IV A, this imaginary charge 
should vanish after functional integration over the time- 
like component of the gauge field Aq (x) . 

The net charge and current are obtained by summing 
the local charge and current densities over fixed time 
slices in Euclidean space: 

J (x ) = Jd 3 x(j (x ,x)) , (28) 

J B (x ) = Jd 3 x(j(x ,x))-B/\B\. (29) 

where 

(jfi (xo,x) ) f = Tr {p- 1 (xo,x;x ,x)7 AI ) . (30) 

We introduce this volume integration as an estimate of 
the macroscopic currents, which are studied in experi- 
ments at scales which are much larger than the typical 
correlation length of the QCD vacuum [2]. 

In Fig. 8 we plot the net overall current Jb (xq) as well 
as the imaginary part of the charge density ImJo(a:o) 
for the 16 3 x 4 lattice, for magnetic fields both parallel 
and perpendicular to the caloron symmetry axis (which 
is parallel to the S3 basis vector) and for different field 
strengths. Despite the fact that we use a fixed and suf- 
ficiently smooth gauge field configuration, the inversion 
of the overlap Dirac operator still takes a lot of com- 
puter time, and the calculation of the Dirac propaga- 
tor D~ l (x, x) in (27) and (26) for all lattice sites would 
be prohibitively expensive. To overcome this difficulty, 
we perform the integration over the spatial coordinates 



x in the definitions of Jo (xq) and Jb (xq) by a Monte- 
Carlo method, spreading the points x uniformly over time 
slices and summing up the local current densities in these 
points. We use from 50 to 100 random points, depending 
on the lattice size. Error bars on Figs. 8 and 9 illustrate 
the estimated uncertainties of the Monte-Carlo integra- 
tion. 

For the parallel case, one can clearly see that the imag- 
inary part of the net electric charge grows with the field 
strength and practically does not depend on xq. For 
the perpendicular case, the charge is close to zero for 
all xq. The net current in the direction of the magnetic 
field Jb (xq) is smaller than the imaginary part of the 
charge density by more than two orders of magnitude. 
The electric current is nonzero in both the parallel and 
the perpendicular cases, but shows an oscillating behav- 
ior, so that the sum over all xq is zero within the errors 
of Monte-Carlo integration. Note that the caloron we use 
has almost static topological charge and action density, 
whereas the gauge fields in its constituent monopoles' 
cores are gauge-rotated in time with respect to each other 
[1 I.]; this twist could cause the change of sign of the 
current Jb (xq). Such behavior of Jb (xq) suggests that 
no net electric charge is transported in the direction of 
the magnetic field in the tunneling event and thus this 
caloron cannot be considered as a model field configura- 
tion for the Chiral Magnetic Effect. On the other hand, 
the induced imaginary electric charge should contribute 
to the fluctuations of the net macroscopic charge, as dis- 
cussed in Subsection IV A above. In the following, we 
consider only the more interesting case when the mag- 
netic field is parallel to the caloron axis (that is, parallel 
to e 3 ). 

In Fig. 9 we plot the imaginary part of the net elec- 
tric charge Im Jo (10), averaged over all xo, on differ- 
ent lattices as a function of the magnetic field strength. 
The value of the magnetic field strength eB on differ- 
ent lattices is kept constant (the total magnetic flux 
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FIG. 9: Imaginary part of the net electric charge Jo (xq), av- 
eraged over all xo, as a function of the magnetic field strength 
on different lattices. 

$b is changed correspondingly). At moderately small 
field strengths the charge rises approximately linearly, 
and then seems to saturate at strong magnetic fields. 
Such behaviour is in complete agreement with the con- 
siderations of the previous Subsection IV A. The solid 
line on Fig. 9 is a linear fit through data points with 
eB < 0.2, which yields J = (8 ± 2) eB. The results ob- 
tained on different lattices agree within the uncertainties 
of the Monte-Carlo integration, which suggests that the 
induced charge is not a finite- volume nor a discretization 
artifact. Some specific lattice artifacts for the density of 
induced charge will be discussed below. 

Next we investigate the longitudinal profile of the 
imaginary part of the charge density Im ( j {x) ) along 
the caloron symmetry axis with x\ = L/2, x-i = L/2, and 
the transverse profile in the direction perpendicular to it. 
We consider the fixed time slice xq = 0, where the net 
charge is maximal. In Fig. 10 we show the longitudinal 
profiles of the charge density on lattices of different size 
and spacing. For the left and the right plots of Fig. 10, 
the magnetic field strengths differ by a factor of two. The 
coordinate X3 on the plots is shifted by ±2 or rescaled by 
2 /3 in order to match the positions of the monopole and 
the anti-monopole on different lattices. These positions 
are marked by red (monopole) and blue (anti-monopole) 
vertical lines on Fig. 10. In what follows, we denote the 
corresponding 2; 3 coordinates as x 3 m ^ and x 3 , respec- 
tively. For the 16 3 x 4 and 24 3 x 6 lattices x ( 3 m) = L/A 

and x 3 a) = 3L/4, For the 20 3 x 4 and 24 3 x 4 lattices we 
simply add more space around the caloron, so that the 
distance between the monopole and the anti-monopole 

Xg°' — Xg" 1 " 1 = L/2 and the position of the midpoint be- 
tween them (x 3 a ^ + x^j /2 = L/2 are fixed. 

One can see that the charge density has a character- 
istic three-peak structure. Namely, there is an excess of 
positive imaginary charge peaked around the monopole 
position, and an excess of negative imaginary charge near 



the anti-monopole. Furthermore, the charge density has 
a pronounced positive peak right in the middle between 
the monopole and the anti-monopole, which is higher 
than the peaks near the monopole and the anti-monopole. 
The height of the middle peak becomes almost two times 
larger as the magnetic field strength doubles. The side 
peaks become also higher and more pronounced. Com- 
paring the charge profiles for all four different lattices, 
we conclude that boundary effects and finite-spacing ar- 
tifacts are rather small and do not exceed 20%. 

In order to investigate the transverse profiles of the 
charge density, we set £3 = x 3 m \ 23 = (xg™* + x 3 a ^ /2 

and x 3 = x 3 and change the transverse coordinate x\. 
The dependence of the charge density on X\ for these 
three values of 23 is illustrated in Fig. 11. The caloron 
symmetry axis goes through the centers of all the plots. 
Again, the coordinates were shifted and rescaled in order 
to match the caloron axis at xi — L/2, X2 — L/2 on 
different lattices. The magnetic field strengths are the 
same as for the left plot of Fig. 10. 

For all three values of X3, the charge density is peaked 
around the caloron symmetry axis. There is also a large 
excess of positive imaginary charge density close to the 
boundaries of the lattice, which clearly diminishes as the 
lattice size increases. A comparison of the profiles for 
the 16 3 x 4 lattice and the rescaled 24 3 x 6 lattice shows 
that the profiles practically do not depend on the lattice 
spacing. This suggests that these boundary peaks are 
finite- volume, rather than discretization, artifacts. For 
asymptotically large volumes one can therefore expect 
that the charge density is localized near the caloron axis 
only. 

These findings can be readily interpreted along the 
lines of the previous Subsection IV A. Namely, the 
monopole and the anti-monopole which constitute the 
caloron create a chromomagnctic field between them. In 
the vicinity of a mid-point this field is approximately con- 
stant, uniform, parallel to the caloron axis and propor- 
tional to the t 3 matrix in the color space. The nontrivial 
holonomy of the caloron solution, i.e. the Polyakov loop 
outside of the monopole and the anti-monopole, is also 
proportional to t 3 and thus plays the role of an imag- 
inary isospin chemical potential. If a constant Abelian 
magnetic field parallel to the caloron symmetry axis is 
superimposed on such a field configuration, an imaginary 
electric charge should appear between the monopole and 
the anti-monopole. 

The origin of the charge excess which is localized di- 
rectly on the monopole or the anti-monopole is not so 
clear to us. Possible explanation could be the follow- 
ing: since the monopole and the anti-monopole also cre- 
ate chromoelectric fields E ~ ±r/r 3 r 3 , in a state of 
thermal equilibrium the fermions would tend to decrease 
this electric field by forming Debye screening clouds 
of negative and positive r 3 -charge, respectively. This 
screening clouds should be seen in the isospin density 
( ip (x) 70I" 3 ?/' {%) ) (note that this is written in the par- 
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FIG. 10: Profiles of the charge density along the symmetry axis of the caloron on lattices of different sizes for fixed values of 
the magnetic field strength eB. Red and blue vertical lines mark the positions of the monopole and the anti-monopole. 
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FIG. 11: Profiles of the imaginary part of the charge density Im ( jo (x) ) in the direction perpendicular to the symmetry axis 
of the caloron for different lattice parameters and at different x%: at the monopole position (left), at the midpoint between the 
monopole and at the anti-monopole (middle) and at the anti-monopole position (right). The magnetic field strength corresponds 
to $s = 5 units of magnetic flux through the 16 3 x 4 lattice. 



ticular gauge where the holonomy is in the redirection) . 
However, they would not contribute to the electric charge 
density as there would be an equal amount of particle as 
well as anti-particle states due to the T3 symmetry. This 
argument, however, is no longer valid when the magnetic 
field is in place, as we have seen that the combined effect 
of magnetic field and nontrivial holonomy is to lift this 
symmetry. This would cause otherwise invisible Debye 
clouds to be seen in the charge density as bumps of pos- 
itive/negative charge around the cores of the monopoles. 

C. An estimate of possible experimental 
importance 

Let us now roughly estimate the importance of the 
described phenomenon for heavy-ion collisions. As the 
simplest approximation, we first adopt the picture of 
"spaghetti" QCD vacuum [31], assuming that it con- 
sists of domains with roughly constant chromomagnetic 
field F. According to (25), the relative decrease of the 
dispersion of charge fluctuations is given by the ratio 

We estimate the typical strength of chromo- 



magnetic field in QCD at temperatures close to the de- 
confinement phase transition as F ~ 0.1 GeV 2 basing on 
the lattice measurements of the finite-temperature gluon 
condensate [32]. As a moderate estimate of the magnetic 
field strength, we take eB = m 2 « 0.02 GeV 2 which 
corresponds to Au-Au collision with impact parameter 
b = 4fm at yfs = 200 GeV [33, 34]. Then from (25) it 
follows that the effect of magnetic field of such strength 
is to suppress the charge fluctuations by several percent. 
This estimate of course does not take into account any in- 
homogeneities of the magnetic and the chromomagnetic 
fields as well as different orientations of the latter in the 
color space and thus can only be considered as a qual- 
itative prediction. Also higher magnetic fields will lead 
to much stronger suppression of charge fluctuations - ac- 
cording to (25), the most significant suppression occurs 
when the strengths of external magnetic field and of chro- 
momagnetic field are nearly equal. The effect might be 
also somewhat larger for real QCD due to larger degen- 
eracy of states for SU (3) gauge group as compared to 
SU (2) gauge group considered here. 

As a more advanced model of QCD vacuum, one can 
consider a model of a dilute caloron gas [15]. We first 
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fix the lattice spacing a in physical units by using the 
expression for the distance between the monopole and 
the anti- monopole in the caloron d = irp 2 T 1 where T is 
the temperature and p is the caloron size parameter, for 
which we take a characteristic value p = 0.33 fm from 
[15]. Taking into account that for the 16 3 x 4, 20 3 x 4 
and 24 3 x 4 lattices the distance between the monopole 
and the anti-monopole is 8 lattice spacings, we obtain 
a pa 0.10 fm for these lattices. Then from the linear fit 
on Fig. 9 we roughly estimate the total charge induced 
on a single caloron as Im J (eB) sa ceB, where eB is the 
projection of the magnetic field on the caloron symmetry 
axis and c pa 2GeV~ 2 . For a dilute gas of calorons and 
anti-calorons with concentration n pa lfm [15] in a four- 
volume V^d we estimate the contribution of the induced 
imaginary charge Im Jo to the expectation value of the 
squared total charge Q as 

6Q 2 ^-l/3((J 2 {eB)))nV id , (31) 

where ((...)) is the average over the caloron ensemble. 
Here we have assumed that the topological charges of 
all "quasi-particles" in a gas are statistically independent 
and take the values ±1 (corresponding to either a caloron 
or an anti-caloron) with equal probabilities. The factor 

7T 

1/3= — / 2irdcos9 cos 2 6> 
4tt J 

— 7T 

comes from averaging over all possible caloron orien- 
tations in the three-dimensional space. For our sim- 
ple estimate, we take V\ d as the product of a typical 
three-dimensional volume of the fireball created in Au- 
Au heavy-ion collisions, V^d ~ (5fm) 3 , times the inverse 
critical temperature of the deconfinement phase transi- 
tion T c w 200 MeV pa (lfm) -1 . The magnetic field is 
again estimated as eB pa 0.02 GeV 2 . With such collision 
parameters, we estimate the contribution of an imaginary 
charge induced on calorons to the fluctuations of the to- 
tal charge as 5Q 2 sa —0.05. It can be easily changed 
by an order of magnitude, for example, if one takes into 
account the non-Poisson distribution of the number of 
calorons in a fixed volume [15]. 

We thus conclude that within the simple models which 
we have considered the effect of the magnetic field would 
be to decrease charge fluctuations in off-central heavy- 
ion collisions by several percent. Since the overall spatial 
volume V^d of the system also decreases for off-central 
collisions, such decrease of charge fluctuations should be 
most clearly seen in centrality dependence of the volume- 
independent observables such as D = ( Q 2 ) / ( N to t ) [35] , 
where Ntot ~ V is the total number of detected charged 
particles. 

We should also make a common cautionary remark 
concerning the use of the grand canonical ensemble ap- 
proximation to describe the fluctuations of charge in 
heavy ion collisions. Namely, the total charge of all pro- 
duced particles is of course conserved, but in experiment 



charged particles are only detected in a certain rapid- 
ity range. The charge of detected particles does fluc- 
tuate and for some reasonable choice of rapidity range 
these fluctuations can be described by a grand canonical 
ensemble of quark-gluon plasma at some finite chemical 
potential. 



V. CONCLUSIONS 

In this paper we have described a mechanism by which 
a combination of magnetic field and isospin magnetic field 
can lead to the emergence of a net electric charge in a gas 
of Dirac fermions in the presence of isospin chemical po- 
tential or a non-trivial holonomy of a non-Abelian gauge 
field. If the system is described as having a constant 
chemical potential, that is, in terms of a grand canoni- 
cal ensemble, the generated charge is borrowed from the 
environment to which the system is coupled and which 
keeps the chemical potential constant. We have also 
demonstrated that in case when no charge is exchanged 
with the environment, the charge density still appears 
due to spatial charge separation, so that the region with 
nonzero field strength acquires the excess of charge of 
one sign, and its exterior acquires the charge of opposite 
sign. In its simplest form the effect can be realized and 
observed in strained monolayer graphene, where a com- 
bination of mechanical strain, magnetic field and valley 
imbalance should lead to the emergence of nonzero elec- 
tric charge. 

A more subtle case is that of non-Abelian gauge the- 
ory, in which the role of the isospin quantum number is 
played by the color of quarks and the nontrivial holonomy 
of gauge field configurations (Polyakov loop) is equivalent 
to an imaginary isospin chemical potential. In this case 
the induced charge is also imaginary for a fixed field con- 
figuration and vanishes after averaging over field config- 
urations with all possible values of the holonomy. How- 
ever, this imaginary charge still yields a negative con- 
tribution to the expectation value of the squared total 
charge in the system. We have explicitly studied the 
induced imaginary charge and its effect on charge fluctu- 
ations in the backgrounds of a constant chromomagnetic 
field and an SU (2) caloron with superimposed Abelian 
magnetic field. We found that the induced charge is local- 
ized between the monopole and the anti-monopole which 
constitute the caloron. Simple estimates made in Subsec- 
tion IV C suggest that the described effect should tend to 
slightly suppress charge fluctuations in off-central heavy- 
ion collisions. 
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